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Abstract 

An action having an acceleration term in addition to the usual ve- 
locity term is analyzed. The quantum mechanical system is directly 
defined for Euclidean time using the path integral. The Euclidean 
Hamiltonian is shown to yield the acceleration Lagrangian and the 
path integral with the correct boundary conditions. Due to the accel- 
eration term, the state space depends on both position and velocity - 
and hence the Euclidean Hamiltonian depends on two degrees of free- 
dom. The Hamiltonian for the acceleration system is non-Hermitian 
and can be mapped to a Hermitian Hamiltonian using a similarity 
transformation; the matrix elements of this unbounded transforma- 
tion is explicitly evaluated. The mapping fails for a critical value of 
the coupling constants. 

1 



1 Introduction 



The action with acceleration arises in many diverse fields; the Euclidean ac- 
tion and path integral, studied by Hawking and Hertog [6j, arises in quantum 
gravity and D-brane dynamics; Bender and Mannheim have extensively stud- 
ied the problem of ghost states in quantum mechanics using the Minkowski 
action [9]; the action appears is the study of mathematical finance, describ- 
ing the time dependence of financial instruments pQ, and the Euclidean path 
integral has been studied extensively in [2], [13]; the action is as an example 
of higher derivative Lagrangians [12j; the path integral for pseudo-Hermitian 
Hamiltonians has been studied starting from the propagator for the theory 
m and HD]. 

Euclidean time is defined by the analytic continuation of Minkowski time 
tM = —iT, where r is Euclidean time. The degree of freedom x does not 
change in going from Minkowski to Euclidean time but Minkowski velocity 
vm picks up an extra factor of i since it is related to Euclidean velocity v by 
the following^ 

doc , doc , / -I \ 

V = -— = I—— = IVm (1) 
ar diM 

The non-Hermitian Euclidean Hamiltonian and path integral are well be- 
haved with a single analytic continuation of time, and no further complexi- 
fication of the degrees of freedom is required. In contrast, an analytic con- 
tinuation of the degrees of freedom is required for the Minkowski case, as 
discussed in [5]. 

^The minus in the definition of Euclidean velocity follows from v = ivM, which is 
consistent with r = Um- 
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2 The Feynman path integral 

Consider the Euclidean time Lagrangian with acceleration given by 
with the acceleration action for finite Euclidean time r given by 



S[x] = / dtC (3) 
Jo 

The Feynman path integral for finite Euclidean time is given by 



lC{xf,Xf-Xi,Xi) = DX e-^l^l (4) 

J {xi,Xi;Xf,Xf) 

DX=AfY[ dx{t) (5) 

*_n J —OO 



where TV" is a normalization constant. The paths have the following four 
boundary conditions 



dx(0) 

x[Oj = Xi] = Xi mitial position and velocity (6) 

x{t) = Xf] — = fiiicil position and velocity (7) 



Since the path integral given Eq. H] is quadratic, it can be evaluated 
exactly using the classical solution. Let Xc{t) be the classical solution given 
by 

6x{t) ^ ' 
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that satisfies the following boundary conditions 



Xc{0) = Xi] — ^ — = Xi initial position and velocity (9) 
Xc{t) = Xf] — = 2;/ final position and velocity 



Consider the following change of integration variables, from x{t) to ^{t) 

x{t) = x,{t) + ^{t) (10) 
with boundary conditions for given by 



^(0) = 0; = initial position and velocity (11) 

^(r) = Xf] — — = final position and velocity 

(JjV 



The change of variables yields 

S[x]=S[x, + ^]=S[x,]+S[^] 

j DX = A/'(r)e'5[^d (12) 

where ^^{T) = J e"^'^! 
Hence, From Eqs. H] and [12] 

IC{xf,xf,Xi,Xi)=Ar{r)e^^'"^'^ (13) 



The evolution kernel given in Eq. [13] has been evaluated explicitly in [6] and 
by solving for the classical solution xdt) and then obtaining S[xc] and 
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Afir). 

As can be directly verified from tlie classical action iS[xc], the classical 
solution Xc(t) given by Eq. [S] yields another equally valid classical solution 
Xc(t) given by the following 

= x,(r - t) (14) 

with boundary conditions, from Eq. [HI given by 

Xc{0) = x^{t) = xf, = -^ = -Xf (15) 

Xc{r) = Xc{0) = Xi\ = = (16) 

^ S[Xc] = S[Xc] 

Hence 

The classical action Sc is given in Eq. I121( the solution is seen to have 
the symmetry given in Eq. [T71 

The evolution kernel, from Eqs. [12] and [TTl has the following symmetry 

lC{xf,Xf; Xi, Xi) = ]C{xi, -Xi] Xf, -Xf) (18) 
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3 Euclidean Hamiltonian and Path Integral 



The derivation of /C(x/, xj; Xj, Xj) was done entirely in terms of the coordinate 
degree of freedom x{t) and made no reference to any other degrees of freedom; 
in particular, the velocity degree of freedom v{t) did not appear in the path 
integral derivation. 

One would like to interpret the evolution kernel /C(x/ the 
probability amplitude for a transition from an initial to a final state vector. 
Such an interpretation of course needs both a state space and a Hamiltonian. 

Based on the boundary conditions given in Eq. [HI it can be seen that the 
state space has to have two independent degrees of freedom, corresponding 
to the two initial conditions given by the initial position x and velocity x. 
Hence, the state space V of non-Hermitian Hamiltonian is taken to have two 
degrees of freedom, namely a position x and a velocity v degree of freedom . 
The Hamiltonian has to be chosen in such a manner that the velocity degree 
of freedom v is constrained to be equal to the velocity x of the coordinate 
degree of freedom x. 

Consider two independent degrees of freedom x and v. The completeness 
equation for the basis states are given by 



A state space representation of the evolution kernel /C(x/, x/; Xj, Xj) is 
derived. It will be shown that the evolution kernel is closely related to the 




(19) 



(x, v\x', v') = 6{x — x')6{v — v') 



6 



probability amplitude of going, in time r, from the initial state \xi,Vi) to the 
final state {xf,Vf \ and is given by 

JCsixf,Vf;Xi,Vi) = {xf,Vf\e~''"\xi,Vi) (20) 

A similar definition is adopted for the path integral of a non-Hermitian Hamil- 
tonian in [8]. 

It remains to be seen as to what is the precise relation of the probability 
amplitude ICs{xf, Vf] Xi, Vi) defined using the state space and Hamiltonian to 
the probability amplitude }C{xf,Xf;Xi,Xi) defined using the path integral. In 
particular, as it stands in Eq. [201 the initial Vi and final velocity Vf have no 
relation with the coordinate degree of freedom x. The Hamiltonian has to 
implement a constraint to set the initial and final state in Eq. |2D]to have the 
same as the boundary conditions given in Eq. |6]for /C(x/, x/; Xj, Xj). 

The Hamiltonian, for infinitesimal time r = e, is given by the Dirac- 
Feynman formula as follows 

{x, v\e-'"\x', v') = c:(e)e^^(^'^''^'^') (21) 

where C(e) is a normalization constant that depends only on e. The discrete 
time Lagrangian is given by 

£(x, x'; V, v') = -I ' - I ' - \[^{x) + $(x')] (22) 

The Minkowski Hamiltonian for the action acceleration has been obtained 
by Bender and Mannheim [5], [9]; they have shown that Hamiltonian and 
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state space for Minkowski time is well behaved, but requires an analytic 
continuation of the degree of freedom. 

The analytic continuation to Euclidean time of the Minkowski Hamilto- 
nian yields a Euclidean Hamiltonian given by 



The term —vd/dx is noteworthy since it does not depend on any coupling 
constant; this term constrains the degrees of freedom and finally leads to the 
constraint that v = —dx/dt, as required by Eq. [H the constraint due to the 
term —vd/dx in the Hamiltonian has a transparent realization in the path 
integral formulation and discussed in Section HI 
Note that H is not Hermitian since 



The transition probability amplitude is given by defining e = t/N and 
inserting — 1 complete set of states given in Eq. [191 Hence, for bound- 
ary conditions given by xq = Xi,vo = Vi]XN = Xf,VN = Vf, the transition 
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amplitude is given by 



]C{xf,Vf]Xi,Vi) = {xf,Vf\e~^^\xi,Vi) 

AT-l N 

= Y\ dXndVnY\{Xn,Vn\e~^"\Xn-l,Vn-l) 
71=1 n=l 

N-1 

= Yi dXn{xN,VN\e~^^\xN-l,VN-l) 
n=l ^ 

'N-1 „ 

JJ / dVn{Xn,Vn\e~^^\Xn-l,Vn-l) 



(25) 



The reason for putting the term {xN,VN\e~''^\xN-i,VN-i) outside the / dvn 
integrations is because, as will be seen in Eq. |271 this term does not depend 
on Vn-1- 

The differential operator H given in Eq. |23l for Xn = x,v = f„; x' = 
Xn-i,v' = Vn-1, yields the following^ 



{x,v\e-'''\x',v') = e-'"^'''^''''^\x\x'){v\v') 

-eH f dp dq j„(x-x' 



ip(x—x')+iq(v—v') 



2-K 2-K 



2n2TT 



C5{x -x + tv) exp{--^(f - v'f - —v'^ - e$(x)} (26) 



Hence Eq. [26] yields the following 



{xn,Vn\e ^■^|x„_i,t;„_i) = C5(x„ - x„,_i + et;„)exp{e£„} (27) 



^Note that the Hamiltonian acts on the dual basis states {x,v\ 
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The 5-function depends only on Vn and results in constraining Vn on reaching 
a boundary of the path integral. 

Eq. [22] yields the remarkable result that, as mentioned earlier, the term 
vd/dx in the Hamiltonian yields a constraint 5{x — x' + ev) on the degree of 
freedom v so that it is constrained to be the velocity, namely v = —dx/dr. 

The appearance of the 5— function in Eq. [27] yields, in the integrand of 
the path integral, the following constraint 



S{Xn - Xn-1 + eVn) ^ Vn = -— (28) 



hm v„ = -^ (29) 
e^o dr 



The Lagrangian £„, from Eqs. [22] and [271 is given by 



lfVn-Vn-l\^ afXn-Xn-l\^ 1 



2V e / 2V e / 2''^^''"^ + '^^''"-'^] 

The path integral and Lagrangian that appears in Eq. [1] makes no ref- 
erence to the integration over the velocity variables. Hence, all the velocity 
integrations J dvn need to be carried out in order that one obtains the ex- 
pression in Eq. [U Remarkably enough, all the j dvn integrations can be done 
exactly due to the 5— function that appears in Eq. [27] 
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Eqs. 125] and 123 yield the following velocity path integral 




n=l 



= C^ / dVnS{xn - Xn-i + eVn) exp{e£„} 




n=l 



(31) 



4 Change of boundary conditions 

The four boundary conditions given in Eq. |9] are solely in terms of the 
position degree of freedom x{t) whereas the boundary conditions given in 
Eq. [201 the defining equation for }Cs{xf,Vf;Xi,Vi) is given in terms of final 
and initial positions and velocities Xf,Vf;Xi,Vi respectively. 

The initial and final time steps in the path integral need to be examined 
carefully to see how the initial and final velocity, vo and vn respectively, can 
be expressed solely in terms of the position degree of freedom. 

Note the integrand of Eq. [311 for n = 1, yields, from Eq. [301 the following 




On performing the J dvi integration, the delta function constrains Vi 
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—{xi — Xo)/e; hence, Ci has the following value 



exp{e>Ci} = exp I - ^ ^a^]_xo _^ _ a_ _ _^ ^^^^| 

= C{e)S{xi — Xq + evo) + 0(e) = C{e)S{xi — Xi + eVi) (33) 

where Xq — Xi and Vq — Vi are the initial position and velocity. 

The final time boundary term for the action yields, using the final velocity 
vn = Vf, the following 

{xn, VNle'^'^lxN-i, vn-i) = C{e)S{xN - xn-i + ^vn) exp{e£Ar} 

= C(e)(5(x/ — xn-i + ef/) cxp{e£Ar} 

since xn — Xf and vn — Vf. 

Collecting all the results yields the discrete time path integral for the 
transition probability expressed solely in terms of the co-ordinated degrees 
of freedom, namely 



lCs{xf,Vf;Xi,Vi) = {xf,Vf\e ''"\xi,Vi) 

N-l „ N 

— C Y\ / dXnS{xf — xn-1 + evf)5{xi — Xi + eVi) exp{e £„} (34) 

n=l n=l 



where C is a normalization. 

The path integral over the velocity degrees of freedom yields, in addition 
to the expected acceleration action, two delta functions. These delta func- 
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tions are crucial in changing the boundary conditions for the path integral 
over the position degrees of freedom. 

The position path integral given in Eq. due to the two delta-functions 
in the integrand, has four boundary conditions for the position degree of 
freedom, namely Xi,Xf due to boundary conditions from the initial and final 
state vectors and two more boundary conditions imposed on xi, x^^i due to 
the two delta- functions resulting from the velocity path integral; in effect, 
these two delta-functions remove two integrations, namely J dxidx^-i in the 
path integral given in Eq. [31] by fixing the value of Xi, x^-i- 

To take the continuum limit define 



x{t) 



dx(t) Xji Xn—l 

dt e 



t = ne 



(35) 



Hence, from Eqs. [33] and [311 respectively 



Xo — Xi 



e 



dx{0) 
dt 



(36) 



6 



dx{T) 
dt 



Xf 



(37) 



From Eq. [31], 




(38) 



Taking the limit of e — )■ using x = [xn — Xn-i]/^ yields 




— - ^{x) 



(39) 
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The delta functions for the boundary values of the functional integral J DX 
are constraints that change the boundary conditions on the path integral, con- 
verting the two position boundary conditions in the path integral J DXDV 
to four boundary conditions for the path integral j DX. Hence, one obtains 
the following continuum result 

K.s{xf,Vf;Xi,Vi) = {xf,Vf\e~''^\xi,Vi) 
DX6{vi + Xi)6{vf + Xf)e^ 

DXe^' 



{xi,Xi=-Vi;Xf,Xf=-Vf) 
= }C{Xi,Xi = -Vi;Xf,Xf = -Vf) 

ICs{xf, Vf] Xi, Vi) = }C{xf, -Vf] Xi, -Vi) 



Recall from Eq. |42l the Hamiltonian for Lagrangian in Eq. |39] is given by 

5 Equivalent Hermitian Hamiltonian Hq 

To analyze H in greater detail, choose the Gaussian potential 

= ^x^ (41) 
that yields the Hamiltonian and Lagrangian given by 



H = -^^-v§^ + \av^ + f ; C = -\ [.x^ + ax^ + ^x^] (42) 
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A parametrization that is more suitable for studying the Hamiltonian and 
state space is given by [5] 

£ = -| [x^ + {cul + ujI)x^ + ululx'^] (43) 

Note the Lagrangian is completely symmetric in parameters cui and UJ2- 
The Hamiltonian given by 

1(9^ ^i7/2| 2\ 2,7222 / A A\ 

if = --^ - + -K + uj,)v + -u.uj.x (44) 

The Hamiltonian acts on a state space V with two degrees of freedom, 
namely position position x and velocity degree of function v. The state vecor 
|\E') is given by 

1^) G V ; {x,v\^) = '^{x,v) (45) 











Complex Branch 




Real Branch 








«^ = Re' 




0)^ = iie"'^ = CO* 




«2 = -Re ' 




-0 2^ 








a 



Figure 1: Parameter branches. 



The parameters have three branches and are shown in Figure [TJ 
• Complex branch a < l^f^. 
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Frequencies ooi , 002 are complex. 

cj^ = uj; = Re'^ : R>0, G [-7r/2, n/2] 

Note G [— 7r/2,7r/2] for all a, /3 > - and this is also the range for 
which the path integral is convergent. 

• Real branch a > 21/^7. 

Frequencies ui,U2 are real and ui > uj2 is chosen without any loss of 
generality. 

ui = Re" ■ U2 = Re'^ : R>0, b e [0, +00] 

For the case of real cui and U2, the entire parameter space is covered 
by choosing say Ui > U2 and the roots are chosen accordingly as given 
below 




uJi > UJ2 for wi, 002 real 



• Equal frequency a = 2a/^. 



uji = 002 ; b = = (j) 



The special case of equal frequency ui = UJ2 is treated in detail in [3]. 
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A similarity transformation is obtained such that 

e-Q/^He"^/^ = Ho (47) 

where the Hamiltonian Hq is a system of two decoupled harmonic oscillators, 
one each for degree of freedom x and v. 

In a pioneering paper, Bender and Mannheim [5] found the operator Q 
for the Minkowski case; the Euclidean version of their result is given by the 
following 

Q = axv- b—— (48) 
oxov 

For the real domain where ui, 002 are real, both coefficients a, b are real and 
hence 

Q = : Hermitian for Ui, U2 real (49) 

The definition of Q continues to hold for the complex domain but Q is no 
longer Hermitian. 

To obtain Hq, the equation for the commutator 

00 ^ 

e-w = ^-[[[o,g],g]...g] (so) 

^ ji— fold commutator 

needs to be applied to O = x, v, d/ dx, d/ dv. 

To obtain the commutator, note that the n-fold commutator of Q with 
X, V, d/dx and d/dv follows a simple pattern that repeats after two commu- 
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tations. In particular, note that 



d_ 
dv ' 



[x,Q]=b 



d_ 



[v,Q]^b 
[-^,Q]^ax, 



[[x, Q],Q] = abx, 



ab—, 
ox 



[[v,Q],Q] = abv, 



ov 



Carrying out the nested commutators to all orders and summing the result 
yields, for a, 6 > 0, the following 



ox 



dv 



/b d 
-sinh(T\fab) — 
a ov 

coshir^/ab)-^ + < sinh(r's/a6)v 
ox \ b 

/b d 
- sinh(r\/a6)— — 

cosh.( rVab)— — h -i/ — smh.(rVab)x 
ov y b 



(51) 



The Euclidean result given above is simpler that the Minkowski commutators, 
which have alternating signs and i's in the various expressions. 

Consider the following equation 



e-Q/^He^/^ = + C,x^ + C,v^ + C,^^ + C,x' + C,v' (52) 



' dv"^ 



dv 



dx dx'^ 



To obtain the factorization of the Hamiltonian into two de-coupled oscil- 
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lators, choose the following values for a and 6, namely 



"-^001002] sinh(vaoj = — 



V^b = \nr^^±^] (53) 



Note the definition of a and b is based on uJi,U2 being real and Ui > U2, and 
which makes the operator Q Hermitian. 
Define 



A = cosh(^— ) 



C 



B 



smh( ) 



W2 



ur — Uo 



- = 'yuiU2 



(54) 



Using the result of Eq. EDand the definitions in Eq. E^lyields the following 
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B 



2, ,2 



1 

1 

'2^ 



cosh 



ab 



sinh 



ab 



(55) 



Similarly, after some simplifications 



C2 



CAB 



2 2 
1^1^ 2 



AB 



7 ^ -\C 

-{A^ + B')+^{ul + ul) 







The constants a and b are chosen so that C2 = C3 = 0; hence one has the 
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following 

C2 = C3 = =^ Determines a and b 



(56) 



The remaining coefficients are given by 



AB J ^ fB 
C4 = — — - + —{uJi + uf 



2 



C, = -^^B^C' + lululA' = lulul 

C, = -ABC + |(a;? + ul)A^ = ^-uj\ (57) 



Collecting all the results yields 



e-«/2^e«/2 = Eo (58) 
^ ^ + + Jc^^^^V (59) 



27^1;^ 2-ibj\dx^ 2 ' 2 

6 The matrix elements of e~^^ 

The Q-operator is given from Eq. BSlby 



92 

Q = axv - b (60) 

oxov 



where, from Eq. 



^ • 1,/ AT\ 2wiW2 

- = 7ix;ia;2 ; smh(VaO) = — ^ 

b ui-^ — 

The finite matrix elements of the e^'^'^-operator are required for many 
calculations involving the state vectors. The exact matrix elements can be 
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obtained by noting that the Hermitian Q-operator factorizes into two decou- 
pled harmonic oscillators by an appropriate change of variables. 
Consider the change of variables given by 



(61) 
(62) 



In these coordinates the a, P sectors completely f actor ize and yield 



1. 1 2 1 1 ^2 

= -2fe + 2^" + 2 V ' 



(63) 



Consider the Hamiltonian of a quantum oscillator given by 



sho 



1 2 2 

2m dz^ 2 



(64) 



with the transition amplitude given by [TT 



K{z,z';t) = {z\e-^^\z') 



271 sinh UT 



exp 



2 sinh(c(;r) 



2; + ) cosh(a;r) — 2zz' 



(65) 



Comparing the a- and /3-sectors of Q with ifsho shows that, for a, b real, the 
a sector is the usual quantum oscillator but the (3 sector yields a divergent 
transition amplitude. The result for the (3 sector is assumed to be given by the 
analytic continuation of the oscillator transition amplitude; this assumption 
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will later be verified by an independent derivation. 

To exploit the quadratic form of the a and /3 sectors, consider extending 
the range of m to the real line. 

a — sector : h = — ; a = mu? =^ a; = ^J~ah (66) 

m 

(3 — sector : b = ; a = —mu'^ =^ co = Vab (67) 

m 

Hence Eqs. |63l |65] and Eq. [67] yield the following 

{a,P\e~^'^\a',P') 

= NaNp exp < — — ^ — [(a^ + a'^) cosh(rv^) — 2a;Q;'] 

I 2osinh(rva6) 

X exp I 4 ^ [(/32 + cosh(rv^) - 2/3/3'] 

I 2o sinh(rv ah) 

= ATir) exp I - ^g{r){a^ + a'^ - - f3'^) + H(r)(aa' - /3/3') 



where 



^(r) = . r coth(rv^) ; ^(r) = . ^ , (6^ 
V V Osinh(rVa6) 



and 



A/-(r) 



a6 / V afe 



27r6sinh(rv a6) y — 27r6sinh(rvct&) 



A/-(r) = r = ^\n{T)\ (69) 

^ ^ 27rV 6|sinh(rv^)| 27r I ^ ^1 ^ ^ 
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Hence a heuristic derivation for Q yields 



{x,v\e-^^\x\v') =M{T)e^-^{-g{T){xv + x'v')+U{t){xv' + vx')} (70) 

The normahzation constant M{t) will be seen to play a crucial role in the 
normalization of all the state vectors of the Hamiltonian H. 

For the real branch, both coefficients Q{t) and 'H(r) are real and hence 
Q = Q'^ is Hermitian. 

Note the form obtained for e~'^^ in Eq. [70] is a major simplification since 
in general, for Hermitian Q, one would need to evaluate 4 + 12 = 16 real co- 
efficients. Instead, the form that has been heuristically derived in Eq. [70] has 
reduced the determination of e~'^^ to that of computing two real coefficient 
functions ^(t) and 'H(r) and a normalization constant A/'(r). 

The operator Q is unbounded and many of the manipulations are only 
formally valid. For example, the identity e~'^^e^^ = I holds for the ma- 
trix elements only in a formal sense. Evaluating the matrix elements of the 
product e~'^^e'^^ yields the following 



{x,v\e-^'^e^'^\x\v') = J d^dC{x,v\e-^'^\^,0{^X\e^^W,v) 
= A/'2(r)e-^(")("'''') j d^dCexp{in{T)i^{v' - v) + in{T)iC{x' - x)} 



27r 



AT^t) ( ] 6{x - - v') = 6{x - x')6{v - v') 
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Hence 



To make the above derivation more rigorous one can analytically continue r 
back to Minkowski time t = —ir, do the computation and then analytically 
continue back to Euclidean time r. This would give the result given above. 

7 Vacuum state; e^*^/^ 

The vacuum state and the dual vacuum state of the Euclidean Hamiltonian 
is obtained using the expression for e*^/^ and e~*^/^ respectively. The vacuum 
can verified to be correct by the direct application of the Hamiltonians H and 

and hence provides an independent verification for the matrix elements 
obtained for Q in Section El 

From the oscillator Hamiltonian Hq given in Eq. [591 the oscillator vacuum 
state, by inspection, is given by 

i/o|0,0) = |(a;i + a;2)|0,0) 

The vacuum state is given by 

l^oo) = 6^5/2 1 0, 0) ; H\^oo) = Eol^oo) ; ^oo = ^[^i + ^2] (72) 
The co-ordinate representation of the vacuum state can be directly ob- 
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tained from the Hamiltonian H and is given by [1] 



^qq{x,v) = (x,t;|^oo) = (a;,t;|e<3/2|0,0) 

= A^ooexp iy--{uji + uj2)uJiUJ2X^ - -(cui + uj2)v^ - ■yuiU2Xvj (73) 

The vacuum state |\l/oo) is real valued and normalizable, with A^qo the nor- 
malization constant. 

It can be directly verified that the dual ground state, from Eq. [73l is 
given by 

^f^o(x, v) = Noo exp + uj2)ujiuJ2x'^ - ^^^i + + l(^i(^2xvj 

More formally, the dual vacuum state obeys the following 

«l = (0,0|e-«/^ 

^g{x,v) = {^^o\^,v) = {0,0\e''^/'\x,v) 
From Eq. [701 r = -1/2 yields 

{x, v\e^^'^\x , v') = 



,^,1, I /" / y/ab\ , . ,, ja(xv' + vx') . 
^^{^) exp < W ^ coth i—-]{xv + x'v') - J } (74) 



sinh , 2 
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Note 



'a , I V cib\ o la 
■ coth 



b \ 2 I \ b gi^i^ ( vJF 



7t^i ; Jv- Y = -fUJi^ujf - uj^ (75) 



^'^'^iV^-r|^^i^"'\^ '''' 

Hence, from Eq. [71] 



Hence 



«(a:,t;|e«/^|e,C)(e,C|0,0) 



where 



•S* = "2 (^iC^ + c^?W2^^ - 2ci^C) - 7C2(Ca; + O + 7Cia;w 



ci = ; C2 = wi y - w| 



Performing the Gaussian integrations over ^ and C yields 
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where 



7 



which is the expected result. 

To determine Nqq note that the Gaussian integration has the following 
matrix 



M = 7 



UJIUJ2 -ci 
-Ci Ui 



\/ det(M) = — UJ2 



Note the determinant of matrix M is negative and hence leads to the correct 
sign for the exponent F of the vacuum state vector. 

Performing the Gaussian integration and using Eq. [76]yields the following 



iVoo = U{ 



1 / 2 3 \ 1/4 

1 / Y^1^2^ 



27r 



27r 



70^1 A / C^I" - UJ2 



VT y v^det(M) 

2 , ,2 / 7V'^2^ 



27r 



i7CJ^''^,/wi - UJ2 



Nqq = {uJiUJ2y^'^i -{Uji + UJ2) 



(77) 



To derive the dual vacuum state, note from Eq. [71] 



{x,v\e '^^'^\x',v') = A/'(-) exp | —'jul^xv + x'v') + 'jui <^ ujI — ' + vx') 



A derivation similar to the one for the vacuum state \E'oo yields the dual 
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vacuum state 



= A^oo exp{--(a;i + uj2)ujiUJ2x'^ - 2 (^i + ^2)v'^ + ^ujiUJ2Xv} (78) 

The norm of a state is defined by the scalar product of state with its dual 
and yields the following norm for the vacuum state 



8 Heisenberg operator equations 

In Schrodinger's formulation of quantum mechanics, all the time dependence 
of a quantum system arises due to the time evolution of the state vector 
I^(^m)), where tM is Minkowski time; the operators O being taken to be 
time- independent. The Schrodinger equation yields the following 




^M) = e-*-^|^) ; {x{tM)\ = {x\e' 



(79) 



with the time-dependent expectation of operator O given by 



E[0;tM] = {x{tM)\0\^{tM)) = (x|e'*-^e»e- 



I*) 



(80) 
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In Heisenberg's formulation of quantum mechanics, all the time depen- 
dence of an operator OuitM) arises due to the time evolution of the operators, 
with the state vector |\E') and the dual state vectors (x| being taken to be 
time-independent; the expectation value is then given by 

E[0-tM] = {x\OH{tMm (81) 
Onit) = e'^^'^Oe-''^'" (82) 

For Euclidean time r, one has the following expressions 



and yields Heisenberg operator equations of motion for Euclidean time 



Or 



The non-Hermitian Hamiltonian H for action with acceleration, from Eq. 
is given by 

jj\ = e-Q/^Hoe'^/^ (83) 
=^H'' = e^'^He^ : pseudo-Hermitian (84) 
^ (e--^) ^ = e--^' = e-«e-^^e'3 (85) 



To obtain the Heisenberg operator equations of motion, it is necessary 
to define the dual state vectors for the pseudo-Hermitian Hamiltonian using 
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Q-conjugation. Namely 

Ix) ^ (x|e-« (86) 

Hence, for Euclidean time, the time dependent expectation value, using 
Eq. |85] yields the following result 

E[0-t] = (x(r)|e-«0|M/(r)) = (x|e-^V-«Oe-^^|M/) 
= (x|e-V^Oe--^|M/) = {x\e-'^OH{rm (87) 

Hence, from Eq. [871 due to the choice of the Hilbert space metric e"*^, 
the Heisenberg time-dependent operator Oh{t) for the pseudo-Hermitian 
Hamiltonian is given by the same expression as for the Hermitian Hamilto- 
nian, namely 

Oh{t) = e^^Ce""^ 
In particular, for the acceleration Hamiltonian given by Eq. [23] 

the Heisenberg equation for the position operator xh{t), with x being the 
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Schrodinger position operator, yields the following 



Mr) = = [H, Oh{t)] = -VH{r) (88) 

Hence, the identification made in the path integral derivation, namely x{t) = 
—v{t), is seen to also hold as an operator equation for the Heisenberg oper- 
ators Xh{t),Vh{t) as shown by Eq. [88l 

To illustrate the role of the Hilbert space metric e^*^ consider the time- 
ordered vacuum expectation value of the Heisenberg operators at two dif- 
ferent (Euclidean) times r > 0. The vacuum state is given by |\l'oo) with 
H'^oo) = i^ool^oo); Eqo = Eq = (wi + uj2)/2] using the rule for forming the 
dual vector of the pseudo-Hermitian Hamiltonian yields the following 

G(r) = (^^o|x^,(r)x^,(0)|^oo) = (^oo|e-«x^,(r)x^,(0)|^oo) 

= ($oo|e"^a;e-"(^-^o)a;|^oo) (89) 

since (v^?ol^^ = ^o(^?ol- 

The propagator G{t) in Eq. [89] is the propagator and is analyzed in detail 
in [3]. 

Recall that the probability amplitude is given in Eq. 

K.s{x^v]x' ,v') = {x,v\e~'^^\x' ,v') 
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and the matrix elements of the the Hilbert space metric e is given (later) 
in Eq. [70] 



(x,f|e '^'^\x' ,v') = N exp {—Q{t){xv + x'v') + 7i{T){xv' + vx')} (90) 

For both the operators e~'^^ and e"'^^, there is no need for an extra metric 
e"'^'^ since Eqs. [20] and [70] are the matrix elements of the operators in a 
complete basis {x,v\ and its dual \x',v'). 

It is only the matrix elements of operators e~'^^ and so on are determined 
for eigenstates and dual eigenstates of H, then the metric is required. 

A verification of the definition of }Cs{x,v]x',v') is given by considering 
the limit of the probability amplitude for r — )■ oo; Eq. [20] yields {Ei is the 
generic energy of the state above the ground state) 

lim )Cs{x,v;x',v') ~ e-^^"" {x,v\e^/^\0,0){0,0\e-^/^\x' ,v') + 0{e-^'^^-^'>'>'>^) 

T—^OO 

= e-^^'"'^ooix,v)^^,ix',v') = e-^^°°^oo(a:,t^)^oo(a:', -v') (91) 

In Appendix B, the value of }Cs{x,v; x' ,v') for large r is obtained from the 
classical solution in Eq. 11271 and confirms the form of the kernel given in Eq. 
[91] is correct and consistent with its path integral definition. 

One can evaluate the matrix elements of e~'^^ in Hilbert space using the 
basis states e*^/^|x,f), the dual basis ((x, f |e'^/^)e"*5 = {x,v\e~'^^'^ and the 
metric e~^, and yields the following 

(x,t;|e-'3/V"^e«/2|^',t;') = (x, t;|e-"^«|x', t;') (92) 
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The symmetry of the matrix elements of (x, f |e~'^''^e~'^^e*^/^|x', f') given 
above in Eq. [92] does not yield the acceleration Lagrangian that results 
for (x, f |e~'^^|x', f'); hence the correct expression for the kernel is given by 
}Cs{x, v; x', v') = {x,v\e~'^^\x' , v'), and which in turn yields the correct matrix 
elements given by the path integral. 

9 Complex cJi, 6^2 

An interesting case for the parameters is the complex domain for which ui, U2 
are complex; the eigenfunctions and eigenenergies are well behaved as long 
as the vacuum state I^E'oo) given in Eq. [73] is normalizable. Let 

uji = Re'^ ; 002 = Re~'^ (93) 
=^ cui + = 2-Rcos(0) ; U1U2 = R'^ (94) 

that yields, from Eq. [721 the vacuum state for the complex domain as 

'^ooix.v) = A^ooexp {-7/?^ cos(0)a;^ - 7i?cos(0)t;^ -^yR'^xv} (95) 

Hence, from Eq. [HSl the vacuum state and Hilbert space for the complex 
domain is well defined for 

cos(0) > =^ -7r/2<0<7r/2 
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10 Many degrees of freedom 

Consider the Hamiltonian that is a quadratic generahzation of the acceler- 
ation action given in Eq. |l2j For degrees of freedom x„, n = 1,2, ...N the 
acceleration Hamiltonian , in matrix notation, is given as follows 



H = -^^S^-S^ -v^ + -vS^'aSv + -xS^pSx (96) 
2dv J dv dx 2 2 ^ ^ ' 



A specific choice is made for the Hamiltonian that is simple and consistent 
with the definitions of v and x given in the defining Eq. [1] The Lagrangian 
is given by 

£ = — - [xS'^'ySx + xS^aSx + xS'^f3Sx'\ (97) 

The (real) orthogonal matrix S and the diagonal matrices are given in matrix 
notation as follow 

SS^ = I ; 7 = diag(7i,72,..,77v) 

a = diag(ai,a2, ••,ajv) ; (3 = diag{(3i, (32, .., (3n) (98) 



Define new variables, in matrix notation 



x = S^z; v = S^u (99) 



The Hamiltonian in Eq. [HS] is given by 



1 ^ 1 92 ^ Q I ^ I ^ 

n=l n=l n=l n=l 
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The corresponding Lagrangian, from Eq. and similar to Eq. |131 is given 

by 

1 ^ 

>C = -- XI ^« [^n + (^In + ^ln)zl + ^L^l^^] (101) 



n=l 

The parametrization is a generahzation of Eq. |16] and is given by 



n 1\FhK\ (102) 

i02n= [\Jan + 2A/7n/3n - \j ttn - '2^J'-fnPn] (103) 

Win > W2n. for Win, real 



The diagonal if given in Eq. |96] is obtained from a Q-operator, which is a 
generalization of Eq. SHI and is given by 

N 

Q = $^Qn (104) 

n=l 

Qn = anZnUn " &n7^ 7{ (105) 

OZnOUn 

with the following values for a„ and 6„ 



/tin / r- , / l^ln + ^2n 

7— = 7nWl„W2n; V «nO„ = In 

On \^ln — ^2n 



(106) 



The diagonal Hamiltonian Hq is given similar to the earlier case. The ground 
state is given by generalizing Eq. [72] and yields the following 

N 

^00(a;, V) = A/"]^ exp |-y (Wi„ + UJ2n)^ln^2nZ'i - ^(Win + UJ2n)u^n " 7nWlnW2n 
n=l 
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In terms of the original coordinates 



"^Ooix, V) ^Afexpl— ^ (^PmnXmXn + QmnVmVn + '^RmnXnVn^ > (107) 
I, m,n=l ) 



where, again in matrix notation 



P = SpS"^] Pn = 7„((^i„ + UJ2n)i^lni^2r. 



11 Conclusions 

The Euchdean acceleration action yields a well defined quantum system that 
has a well defined Euclidean Hamiltonian. The Euclidean path integral is 
the appropriate formulation for numerical and Monte Carlo simulation of 
the system. The path integral using the action and the Hamiltonian yield 
the same result due to a constraint term in the Hamiltonian. 

The state space of the pseudo-Hermitian Euchdean Hamiltonian has a 
state space metric e"*^, which is a natural generalization of the state space 
of quantum mechanics. The metric e"*^ is an unbounded operator that maps 
the non-Hermitian Hamiltonian to the oscillator Hamiltonian; the matrix 
elements of e"*^ were evaluated exactly and shown to yield the expected 
similarity transformation. 

As was seen in various derivations carried out, the results from Minkowski 
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time serve as a useful guideline for the derivations, but given a plethora of i 
and various ± signs that differ between the Euclidean and Minkowski results, 
all the derivations for the Euclidean have to done independently from the 
Minkowski case, and from first principles. 
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A e~^'^ and similarity transformations 

The heuristic derivation for Q{t) and 'H(r) was obtained by an analogy with 
the oscillator Hamiltonian and cannot be assumed to be correct since the /3- 
sector yields an un-stable Hamiltonian. The result needs to be independently 
verify. 

The fundamental similarity transformations by e^'^^ of operators x,d/ dx, v 
and d/dv are directly obtained using the result given in Eq. [70] and shown 
to be identical to the defining equations for Q given in Eq. |511 

Recall from Eq. [5T]that e^"^*^ yields the following similarity transforma- 
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tions 



rv ab) 



d_ 

dv 



I. e '^^xe^^ = coshi T^/ab)x + \ —smh( 

V a 

II. e"^*^— e^*^ = cosh.( Ty/ab)— — \- \ - sinh(r\/a6)f 

ox ox \ b 

/b d 
- sinh(rA/a6) — 

IV. e~'^^-^e'^^ = coshirVab)-^ + W ^ sinh(rA/a6)x 
ov ov \ b 



Consider the operator equation 



e ^'^xe^^ = cosh.{T\/ab)x 



smhiry/ab)-^ 
a ov 



The matrix element of e '^^xe'^^ , using / from above, is given by 



{x,v\e '^^xe^^\x' ,v') = cos\i{t\/ ab)x5{x — x')5{v — v') 



+ \ - smh(TVab)-^6(x - x')5(v - v') (109) 
\ a Ov 



Eq. [70] yields 



{x, v\e '^^\x', v') = N exp{—g{xv + x'v') + h{xv' + vx')} 



(110) 



where 



9 = Q{r) = 
N = ATir) 



a cosh(rv ab) 

b sinh(rv^) 

i 



h = nir) 



b sinh (rv^ofe) 



27r?/(r) 
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The left hand side of Eq. I1U9I yields 



(.,.|e-.e-|.,.)^y.«(..He-|«>«,CI«e-|.'../) 



But 



^exp{-g{xv + ^C) + h{xC + <)} 



— sinfrAi/afo)^ + cosh(rv^a6)x 
a ov 



(111) 



Hence, using e '^^e^^ = I given in Eq. [71], we obtain the expected result 
that 



{x,v\e '^^xe^^\x' ,v') = cosh(rva6) + w ^ sinh(r\/a6)^ {x,v\e ''^e'^^\x' ,v') 



I. e '^^xe^^ = cosh{T\/ab) + \ sinh(r-\/a6) 



A. 

' dv 



Consider 



(x,t;|e--%e-«|x',t;') = / « {x,v\e-^'^\^,0{^X\Ce^''\x',v') 



d^dCC{x,v\e-^''\^X){^X\en^',v') 



b d 

sin]iiT\fab)— — h cosh(rV^)f I 6{x — x )5{v — v') (112) 
a ox 
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and yields the expected result that 



- sinh(rA/a6) — 



Consider the following matrix element 



N / d^dC 



d_ 



(e,C|e^«|x',^') 



d^dC gC-hv {x,v\e-^^\^,0{CX\e^'^\x,v) 



Using Eq. 11121 to replace C in above expression yields, 

e"'^'^— e^*^ = g ^]J^ smh.{T\/ab)— — h cosh{T\/ab)v^ — hv 

d d [h 

=^ III. e"^*^— e^*^ = cosh(rya6)7— + \ - smh(T\/ab)v 
ox ox \ a 



And lastly, similar to above derivation 



{x,v\e-^^-^e^^\x',v') = N f d^dC 
dv J 



(e,C|e^«|x',T;') 



« 9^-hx (x,t;|e--«|e,C)(e,C|e^^|x',t;') 



Using Eq. Illll to replace ^ in equation above yields 



e"^^— e^*^ ^ [ V ~ ^^^^■^(^^'^^)^ + cosh(rva6)x ) — hx 

d d Pb 

IV. e^^^—e"^ = cosh.{T y/ab)— + \ - cosh{T\/ab)x 
ov ov \ a 
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The finite matrix elements of e^^'^ produce all tlie defining similarity 
transformations on tlie x, v degrees of freedom and verifies tliat tlie lieuristic 
derivation of tlie matrix elements of e"^*^ is, in fact, correct. 



B The classical solution 

Consider the following parametrization of the acceleration Lagrangian 



The parametrization chosen in Eq. 11131 is more suitable for studying the 
classical solutions and is different from the one given in Eq. |13i 
The Euler-Lagrangian equation 



yields the equation of motion; the classical solution xdt) satisfies following 
equation of motion 




(113) 



cl'x cit) — 2hxc{t) + cxcit) = 0. 



(115) 



Choose the boundary conditions 



Initial values 



a;(0) 



Xf = xi, x{0) 



Xf = X2 



-Vf 



Final values 




(116) 
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From the definition of tlie probability amplitude given in Eq. 



ICs{xf,Vf;Xi,Vi) = {xf,Vf\e ^ \xi,Vi) = / DXe 



witli boundary conditions given by Eq. 11161 
Define tlie parameters r and co by 



r + iuj 



b + %\J ac — 



(117) 



From Eq. IllSl the classical solution of equations of motion is given by [13 



Xcif) = {di sin ut + 02 cos ut) + e ^ {a^ sin ut + 04 cos ut) (118) 

The parameters ai, .., 04 are obtained from the boundary conditions and given 
by the following. 



ai = r 



r^x/e^'"^ sin(2rci;) + uvje^'^^ — rvje'^^'^ sin(2ra;) + rcoxfe^^^ cos(2ra;) — rcuxf 

— uvf — 2r'^Xie^^ sin(ra;) — 2rt>je^'^ sin(ra;) — we^"^ (e^*"^ — l) cos(rci;) (f j + rxi 

— cj^Xjc'"^ sin(ra;) + u'^XiC^^'^ sin(rci;) 



02 = r 



'Xf\—e )+rvfe + re cos[2tuj) [rxj — vj) — u xje 



2rT 



— rcux/e^'"^ sin(2ra;) + ou'^xj — cuf je*"^ sin(ra;) + uvie^^^ sin^ruj 



I 2 rT { Itt 
+ UJ Xit ' " 



^g2rr _ cos(ra;) + rujXie^'^ sin(rci;) + rujXie''^'^ sin(™) 
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as = Te'' 



^x/e''"^ sin(2ra;) + uvfc"^^ — uvje^"^^ + rvje^"^ sin(2ra;) 



+ rooxje ^'^ — rujXfe^'^ cos(2™) — 2r x^e ^"^ sin(™) — 2rvie '"^ sin(™) 
— u (^e^^^ — l) cos(™) {rxi — Vi) — co'^Xie^^^ sin(™) + tu'^Xi sin(ru;) 



rr 
e 



r'^Xf (— e*"^) — rvfC^^ + re^^ cos{2tuj) {rxj + vj) — co'^Xfe^^ 



+ u'^XfC^^^ + rcoxfc''^ sin(2™) — covie^'^^ sin(™) — co'^Xi (e^^"^ — l) cos(™) 



rojXie ^"^ sin(rci;) — rooxi sin(rci;) + oovi sin(rw) 



In the above equations, F is 



uj2 + ^2g4rr + 2r^e'^'^^ cos(2™) - 2e^'-^ {r^ + u^) 



(119) 



Choosing boundary condition as Eq lll6l the classical action yields 



1 4 

Sc = Sc{Xf,Vf,Xi,Vi) = -7^ ^ XjMijXj 



i,j=i 



From Eq. [T71 the classical action has the symmetry 



Sc[Xf, Vf] Xi, Vi] = S[Xi, -Vi]Xf, -Vf] 



(120) 



Hence under this boundary condition, the solution of matrix M satisfies the 
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following symmetry 



Mn = M44, M22 = M33 

M12 = -M34, Mi3 = -M24 



and the action can be simplified to 



1 " ' 2 , „2\ 1 71^ /„.2 , „.2> 



t;/, t;^, Xi) = --Mn{xf + xj) - -M22K + ^;^) 
+Mi2{xiVi - XfVf) + Mi3(a;if/ - XfVi) - MuXiXf - M2zViVf (121) 



Solution of Mjj is listed below. 



Mil = r 2aruj (r^ + w^) (g^rr _ 1) + 2re^''^ sin(2™)) 
M12 = r L^e^''^ (2ar2 - 6) - 2r2e2'-^ (200;^ + h) cos{2tuj) 

-u" {b-2ar^)+2be^''^ {r^ + u^) 
Mi3 = r [4ara;e^^ (e^^^ - l) (r^ + u'^) sin(Tu) 

Mi4 = r [ - Aarue'^ (r^ + uo^) (r (e^"" + l) sin(rc<j) + 00 {e^'^ - 1) cos(™)) 
M22 = r [ - 2aruj {u {-e^'^) + 2re^''^ sin(2ra;) + u] 
M23 = r Uaruje"^ (r (e^"^ + l) sin(ra;) - uj {e^^^ - l) cos(ra;)) 
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B.l Infinite time limit of classical action 

The infinite r limit of matrix elements M under the boundary condition 
Xf,Vf,Xi,Vi are given by the following. 



lim Mil = 2ra(r^ + w^) ; lim M22 = 2ra ; lim M12 = 2r^a - b (122) 



lim Mi3 = lim Mm = lim M23 = 



From the definition of the probability amplitude given in Eq. [20] 

}Cs{xf,Vf;Xi,Vi) = {xf,Vf\e~'^^\xi,Vi) = {xf,Vf\e^f'^e~''^°e~^f'^\xuVi) 
yields, for infinite r limit as given in Eq. [91], the following 

lim Ks{xf,vf,Xi,Vi)-e-'^°'{x},Vf\^oo){'^oo\xuV.:) + 0{e-^^''^-'''^^) 

T— > + 00 

Hence, from Eq. I123[ the infinite limit of transition amplitude is the product 
of vacuum state \E'^(xi,fi) and \E'oo(a;/, f/) and the classical solution yields 

lim K{xf,Vf,Xi,Vi) = lim AT'e^^^^^''^^'^-^'^ 

r— ^+00 T— >+oo 

= e-^^o°^oo(a;/, Vf)^oo{xi, -Vi) (123) 



The infinite time classical action has precisely the factorization of the 
initial and final boundary values as given in Eq. 11231 Therefore, the vacuum 
state, obtained from Eqs. 11211 and 11221 and, defining xj = x,Vf = v, yields 



45 



the following 



*oo(a;,^^) = lim A/'exp ( -^Mnx^ - ^Maat;^ - M12XV 

T-5> + 00 V ^ ^ 

= A/'exp^— ra^(r^ + u'^)x^ — rav^ — (2r^a — h)xv^ (124) 



where M is fixed by normalizing VL{x, v). 

To make connection with the earlier parametrization, rewrite the La- 
grangian in uji and UJ2 parametrization (the only difference with the parametriza- 
tion given in Eq. |43] is in the overall factor of a instead of 7) 



C = — a + 2-x^ + -X' 
2 \ a a 

= -^a(x^ + (wi + u;^)i'^ + (^1^2)2;^ ) (125) 



where ui and U2 arqf 



wi = ( \/ b + y/ac + \/ b — y/ac 



C02 = I V + \/ ac — \i — yj ac 



The definition variables r and u in terms of ui and is given by 



uoi = r + iuj ; uj2 = r — iuj 



3 



Note that ui and tj2 in Eq. 1461 have a different parametrization, given by 
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The two parameterizations have the following relationship 

+ ^2 _ ^2^2 . ^ = 1 + ; b = f K^i) (126) 

Eqs. 11241 and 11261 hence yield - replacing a by 7 to conform to the notation 
given in Eq. |l6]- the vacuum state given in Eq. [731 namely 

^oo(a;, v) = A/'exp y--{uji + uj2)ojiuj2X^ - ^ (^1 + ^2)^ ^ - 'yujiijj2xvj (127) 

The definition of the evolution kernel lCs{xf, Vf; Xi, Vi) = {xj, Vf\e~'^^\xi, Vi) 
given in Eq. [20] is seen to be correct since the evolution kernel obtained from 
the classical solution also gives the same result for the vacuum state as given 
by the Hamiltonian in Eq. [731 
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